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ABSTRACT 

We present new Hubble Space Telescope images of high-velocity H-a and 
Lyman-a emission in the onter debris of SN 1987A. The H-a images are domi- 
nated by emission from hydrogen atoms crossing the reverse shock. For the hrst 
time we observe emission from the reverse shock snrface well above and below 
the eqnatorial ring, snggesting a bipolar or conical strnctnre perpendicular to the 
ring plane. Using the Ha imaging, we measure the mass flux of hydrogen atoms 
crossing the reverse shock front, in the velocity intervals (—7,500 < Vobs < —2,800 
km s“^) and (1,000 < Vobs < 7,500 km s“^), Mh = 1-2 x 10“® Mq yr“^. We also 
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present the first Lyman-a imaging of the whole remnant and new Chandra X- 
ray observations. Comparing the spatial distribution of the Lyman-a and X-ray 
emission, we observe that the majority of the high-velocity Lyman-a emission 
originates interior to the equatorial ring. The observed Lyman-a/H-a photon 
ratio, (R{La/Ha)) ~ 17, is signihcantly higher than the theoretically predicted 
ratio of ~ 5 for neutral atoms crossing the reverse shock front. We attribute 
this excess to Lyman-a emission produced by X-ray heating of the outer debris. 

The spatial orientation of the Lyman-a and X-ray emission suggests that X-ray 
heating of the outer debris is the dominant Lyman-a production mechanism in 
SN 1987A at this phase in its evolution. 

Subject headings: supernovae: individual (SN 1987A) — shock waves — circum- 
stellar matter 

1. Introduction 

The reverse shock (RS) in SN 1987A is the surface where the freely expanding debris 
is suddenly decelerated as it encounters circumstellar matter. Hydrogen atoms crossing the 
shock are excited by collisions in the shocked plasma and emit H-a and Lyman-a photons 
having Doppler shifts corresponding to the projected velocity of the freely expanding debris 
immediately inside the shock. This emission was predicted by Borkowski et al. (1997) and 
subsequently observed and interpreted by Sonneborn et al. (1998) and Michael et al. (1998). 

The freely expanding supernova debris obeys a “Hubble’s Law” , i.e., the Doppler velocity 
is given by Vobs = z/t , where .2 is the projected distance along the line of sight measured 
relative to the center of the explosion and t is the time since the supernova explosion (f« 27.5 
years in mid-2014). Thus, surfaces of constant Doppler shift are planar sections of the debris. 
The fluxes of H-a and Lyman-a photons emitted at the RS are directly proportional to the 
fluence of hydrogen atoms across the shock. 

These facts create a unique opportunity to map the three-dimensional structure of the 
RS in SN 1987A and provide an empirical framework for studies of the shock interaction. 
This opportunity has been exploited in observations with the Space Telescope Imaging Spec- 
trograph (STIS). In the STIS observations, the H-a emission from the RS was clearly visible 
near the circumstellar equatorial ring (ER), where the shock is strongest (Michael et al. 
1998, 2003; Heng et al. 2006; France et al. 2010). The ring is inclined ~ 43°, north being 
the near side (Panagia et al. 1991; Plait et al. 1995; Sugerman et al. 2005). Michael et 
al. (2003) modeled the RS surface, concluding that the observed emission was conhned to 


3 


within ~ ± 30° of the ring plane, similar to the opening angle of the ring observed at radio 
wavelengths (Ng et ah 2013). However, the STIS exposures were not deep enough to detect 
H-a emission at high latitudes (i.e., out of the ER plane), where the fluence of atoms is much 
smaller than near the equator. Moreover, none of the spectroscopic observations included 
enough STIS slit locations to map the entire shock surface. 

In this Letter, we describe the results of an observing campaign to map the high-velocity 
H-a and Lyman-a emission using filters that transmit the blue- and red-shifted emission from 
the debris while suppressing the bright emission from the ER. In Sections 2 and 3, the new 
HST observations are described; for the first time we detect emission from the shock above 
and below the ring plane. In Section 4, we compare H-a, Lyman-a, and X-ray observations, 
arguing that the majority of the Lyman-a flux is driven by X-ray heating of the outer ejecta, 
as proposed in our earlier work (e.g., Larsson et ah 2011; France et ah 2011; Fransson et ah 
2013). 


2. HST Observations and Image Reduction 

We obtained images of SN 1987A with the WFC3 camera in four hlters to isolate the 
high-velocity H-a: blue high-velocity H-o (“H-a Blue”, F645N), red high-velocity H-o (“H-a 
Red”, F665N), and two narrow-band hlters centered on shocked emission from the ER and 
hotspots (narrow H-a, F656N and [N II] A 6583 A, F658N). The WFC3 images each have 
exposure times of 2880 seconds and were acquired on 21 June 2014. Figure 1 shows the band 
passes of these hlters superposed on our 2014 STIS spectrum of SN 1987A (see also Larsson 
et al. - in prep.). 

The “H -a Blue” and “H-a Red” images include a contribution from emission lines 
originating on the ER (see Fransson et al. - in prep, for a discussion of the ER emission). 
The narrow band images are used to subtract these contributions to the high velocity material 
images (Figure 2). The F656N image provided the best template for the north-dominated 
high-velocity “H-a Blue” image. The F658N hlter contains a signihcant contribution from 
redshifted H-a on the southern side of the remnant, therefore this hlter provided the best 
template for the south-dominated high-velocity “H-a Red” image. In Figure 2 (right), the 
ER appears black; however this is an artifact of the intentional oversubtraction and stretch 
of the images for display. For the quantitative analysis in Section 4, a smaller ER subtraction 
is applied; the emission is nearly continuous across the north and south ER region. 

We observed the Lyman-a emission from SN 1987A with the ACS/SBC F122M hlter 
for three orbits (9200 seconds; Figures 3 and 4). The F122M hlter bandpass is shown on an 
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HST-COS spectrum (France et al. 2011) in Figure 3. In this case, resonant scattering by 
the neutral ISM of the Milky Way and Large Magellanic Cloud acts as a natural blocking 
hlter for low-velocity emission from the ER and no additional subtraction is necessary. A 
UV-bright B star (“Star 3”) appears to the southeast in all imaging bands, providing a 
common reference point. We complement the HST imaging with a recent (2014 March 19, 
day 9885) 0.3 - 8.0 keV image from Chandra (Frank et al. - in prep.). 


3. High-velocity H-a Imaging 

The “H-a Blue” hlter samples H-a emission with Doppler shifts —7500 < Vots < — 
2800 km s“\ while the “H-a Red” hlter samples -1-1000 < Vobs < +7500 km s“^ (Figure 
2). Regions of interest are labeled A, B, C, and D; “-b” and “-r” denote blue and red H-a 
velocities, respectively. Emission from the A-b and A-r regions are produced by H atoms 
passing through the RS near the north and south side of the ER, respectively. The bright H- 
a emission features (A-b and A-r) from the RS have been partially mapped in 3 dimensions 
by STIS observations (e.g., Michael et al. 2003; Heng et al. 2006; France et al. 2010). 

In the “H -a Blue” image (Figure 2, top) a bright crescent-shaped feature covers the north 
side of the ER (as dehned by the [N II] contours), labeled B-b. Given the 43° inclination of 
the ring, the B-b region must originate from signihcantly above the ring plane for it to be 
visible in this conhguration. This image also shows faint streaks that cross the ring extending 
towards the south (C-b). The features labeled A-r, B-r, and C-r are analogous to the blue 
side, but reflected about the major axis of the ER. The additional feature labeled D-r is 
probably emission from the internal debris, which is expanding with radial velocities ranging 
up to ~ 3500 km s“^ (McCray 1993; Larsson et al. 2013). 

Combining the “H-a Blue” and “H-a Red” images, we measure the spatial extent of 
the high-velocity H-a regions (positional uncertainties < 0.05"). The high-velocity material 
interior to the ring (A regions) has north-south and east-west diameters of 0.94" x 1.31", 
which when including the inclination translates to physical dimensions of (9.9 x 9.8) x 10^^ 
cm (assuming d = 50 kpc). This suggests a symmetric, radial flow with a maximum radial 
velocity of 5.7 x 10^ km s“^ (at 27.5 years since the supernova explosion). The B regions 
on the blue and red side have maximum angular distances from the center of the remnant 
of 0.82" and 0.98", respectively. These correspond to maximum de-projected distances and 
radial velocities of 8.6 x 10^^ cm and 10.0 x 10^ km s“^ for the north/blue and 10.4 x 10^^ 
cm and 11.9 x 10^ km s“^ for the south/red. The C region has a maximum observed 
angular extent of 1.37" from the center of the remnant, corresponding to 14.5 x 10^^ cm and 
16.7 X 10^ km s“^ assuming free expansion. 
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3.1. Morphology of the Reverse Shock Surface 

Fransson et al. (2013) noted that the large observed velocity of the H-a at late times 
{Vobs ^ 11 X 10^ km s“^ at t > 20 years) suggests that the reverse shock front has expanded 
beyond the ring plane. They conclude that the RS extent in the polar direction is ~ 15% 
larger than the ER radius and ~ 40% larger than the RS radius in the ring plane. We hnd 
de-projected velocities of 10 - 12 x 10^ km s“^ in the B-b and B-r regions, corresponding 
to RS radii 40 - 70% larger than the 6.1 x 10^^ cm ring radius. The C regions extend to 
2 - 2.5 times the equatorial ring radius. 

To interpret the high-velocity H-a images, we recall that surfaces of constant Doppler 
shift are planar sections of the freely expanding supernova debris. Therefore, the pass bands 
of the “H -a Blue” and “H-a Red” hlters map to slabs in physical space, as illustrated in 
Figure 5. Portions of the RS surface visible through the “H-a Blue” hlter (colored blue) 
reside within the slab delineated by the blue vertical lines, while those visible through the 
“H-a Red” hlter (colored red) reside within the slab delineated by the red vertical lines. The 
density of the outer debris crossing the RS front, and therefore the mass hux, is psn oc po 
t~^ V~^ (Luo et al. 1994). The much higher density of the debris crossing the RS in the 
equatorial direction relative to the polar direction {per/ P poi ~ 10“^) makes the A region much 
brighter (detectable spectroscopically within ~ 10 years of the explosion). The lower density 
material howing out of the ring plane (B and C regions) is fainter, only now visible with the 
dedicated high-velocity H-a imaging and sufficient time for the polar shock structure to be 
angularly separated from the ring plane. 

While much of the RS surface is not imaged because of the discrete velocity surfaces 
transmitted through the hlters (shaded gray in Figure 5), the emerging picture for the ex- 
tended H-a emission in the off-plane regions suggests a bipolar morphology. This is expected 
due to the high pressure of the equatorial plane, owing to the presence of the ring (Blondin 
et al. 1996). Assuming the outer ejecta pressure is proportional to psn V^, the C polar re- 
gion pressure is ~ 1/1800 of the equatorial region pressure. This general picture resembles 
the bipolar surfaces suggested for the pre-supernova interacting wind picture by Blondin & 
Lundqvist (1993), however, recent studies of a ‘pre-explosion twin’ to the SN 1987A system 
(SBWl; Smith et al. 2013) suggest that interacting winds may not be required to explain 
the nebular morphology. 

A second possible explanation for the C region emission is a conical or helical structure 
(e.g. Smith et al. 2007, 2013) connecting the ER with the extended outer rings (ORs). 
In this scenario, 1 ) high-velocity ejecta crossing this structure would be limb-brightened, 
explaining the bright edges in the C-b and C-r emission, and 2) the southern C-b streaks 
would be offset to the east compared to the northern C-r streaks by an amount comparable 
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to the E-W offset of the ORs (e.g, Sugerman et ah 2005), which is indeed what is observed in 
Fignre 2. In this picture, the C emission would essentially be a new SN 1987A intermediate 
ring system. Deep H-« imaging spectroscopy to map the region could distinguish between 
these scenarios and may shed light on the progenitor structure, including the possibility of 
a binary merger prior to the supernova explosion (Morris & Podsiadlowski 2007). 


4. High-velocity Lyman-a Imaging 

Figure 3 illustrates the contribution of different far-UV emission components to the 
F122M image. Interstellar H I blocks the center of the Lyman-a line prohle over the velocity 
range —1500 km s“^ < Vobs < +1500 km s“^ (France et al. 2011), and so excludes emission 
from the shocked ring and hotspots, which have Doppler shifts AV< 300 km s“^ (Pun et al 
2002). Based on the HST-COS spectrum, we estimate that 90 - 95% of the total emission 
in the F122M image is contributed by high-velocity Lyman-a. The remaining 5 - 10% is 
mostly N V A 1240 A emission from the ER and high-velocity ions crossing the RS front. 

Separation of the peak Lyman-a and ER emission is complicated on the bright northern 
side by the 43° projection on the plane of the sky, however, there is separation between the 
Lyman-a and the H-a ring prohles in the East- West direction. The East- West ring diameters 
are 1.44" (± 0.05"), 1.62" (± 0.05"), and 1.63" (± 0.10") for Lyman-a, H-a ER (F656N), 
and X-ray images, respectively. The Lyman-a image is slightly offset with respect to the 
ER, with the largest separation (1.35 x 10^^ cm) at the western edge (Figure 4, bottom). 
While the X-ray resolution is not as high as the HST imaging of the ER, there is a clear 
offset between the X-ray and Lyman-a prohle peaks. Unlike the X-ray and Lyman-a images, 
which are brightest in the NW quadrant, the radio image (Zanardo et al. 2014) is brightest 
in the East. 


4.1. The Lyman-a/H-O! Ratio: X-ray Heating 

Neutral Hydrogen Mass Flux and the Lyman-a/H-a Ratio - If all of the high-velocity 
hydrogen emission originates from collisional excitation as neutral hydrogen atoms cross the 
RS front, the total hux is proportional to the huence of hydrogen atoms across the shock 
boundary. Integrating the individual band images, we calculate the total number of atoms 
in a given velocity range. We take the total (Galactic + LMC) reddening towards SN 1987A 
as a standard ISM curve with Ry = 3.1 and E[B — V) = 0.19 (see France et al. 2011 for 


7 


a discussion). The total rate of hydrogen atoms crossing the RS front is approximately: 

Nh = Ainf X F(AA)^AA x P(H)~^ x Cjsm (1) 

he 

where d ~ 50 kpc (~ 10% uncertainty; Panagia et ah 1991), F(AA) is the image-integrated 
flux density in a given band (~ 5% uncertainty), A A is the effective bandpass of the hlter, 
P{H) is the number of line photons emitted per atom, and Cism is the interstellar attenu- 
ation correction factor. P{H) ^ 0.2 for H-a and 1 for Lyman-a (Michael et ah 2003; Heng 
& McCray 2007). 

For the “H-cc Blue” image, AA ~ 100 A and Cjsm = 1-57 ± 0.08, giving a total hydrogen 
fluence of 4.1 (± 0.7) x 10"^® atoms s“h For the “H-a Red” image, AA 120 A and Cjsm = 
1.53 ± 0.08, giving a total hydrogen fluence of 3.3 (± 0.5) x 10^® atoms s“^. Taken together, 
we estimate the total mass flux of high-velocity hydrogen atoms crossing the RS is Mh = 

7.4 (± 1.2) X 10^2 g s-^ (1.2 X 10-3 M0 yr-^). 

The analogous calculation using the F122M Lyman-a image has AA ~ 100 A and Cjsm 
= 6.7 X 1.63, where 6.7 (± 1.7) is the correction for dust attenuation and 1.63 (± 0.08) 
is the empirically determined correction for the interstellar Lyman-a line core attenuation. 
This calculation shows a Lyman-a-to-H-a photon ratio, R{La/Ha), of ~ 17 ± 6, a factor of 

3.5 larger than attributable to excitation in the RS alone. Therefore, an additional Lyman- 
a emission mechanism is present. Note that the reddening correction does not account 
for differential extinction within the supernova debris, so the intrinsic R{La/Ha) may be 
greater. Lyman-a excess was observed by previous authors using spectroscopic observations 
of small spatial regions of the SN 1987A inner ring region (Heng et ah 2006; France et ah 
2010 ). 

Combining the “H-a Blue” and “H-a Red” images with the Lyman-a map, we hnd the 
baseline R{La/Ha) interior to the ER is 8 - 10, at the center and southeastern sides of 
the inner remnant. R{La/Ha) reaches a maximum on the western and northwestern side 
of the inner remnant where the Chandra emission is brightest (Figure 3). The maximum 
R{La/Ha) is ~ 35 near PA = 270°, just interior to the brightest X-ray emission from the 
ER. R{La/Ha) values ~ 20 are found towards PA ~ 60°, the location of the northeast X-ray 
maximum. 

The Role of X-rays - The observed offset between the X-ray and Lyman-a emitting 
regions is consistent with the expected separation between the ER and the outer extent of the 
RS. Following the observation that the dominant energy input in the unshocked supernova 
debris has evolved from internal heating by radioactive decay to external heating from X-rays 
produced at the ER (Larsson et ah 2011), Fransson et ah (2013) presented calculations of 
the X-ray energy deposition into the ejecta. The heating of the outer debris is dominated by 


X-rays with E < 0.5 keV while higher energy X-rays penetrate deeper to heat the interior. 
X-rays propagating into the outer debris produce fast photoelectrons, which deposit most 
of their kinetic energy as heat, provided that the debris has fractional ionization He/nn 
> 0.03 (Xu & McCray 1991). Radiative cooling by thermal excitation of Lyman-a limits 
the debris temperature to T < 10^ K, at which temperature thermal excitation of H-a is 
negligible^. Therefore, we expect the Lyman-a/H-a ratio to be at a maximum nearest the 
source of the X-ray heating. The spatial stratihcation of the Lyman-a and X-ray emission 
and the enhanced Lyman-a /H-a ratio support the scenario where X-ray heating of the outer 
debris is responsible for the majority of the observed Lyman-a emission from SN 1987A at 
present. 

The total Lyman-a luminosity in the F122M band is = 2.5 (± 0.7) x 10^® erg s“^, 
~ 70% of which is in excess of the expected RS emission. The total flux of the X-ray ring 
is 1.04 X 10“^^ erg cm“^ s“\ or, Lx(0.3 - 8.0 keV) = 3.1 (± 0.3) x 10^® erg s“^. However, 
most of the heating radiation is in the soft X-ray/EUV band (0.01 - 0.5 keV) which is not 
observed at Earth because of interstellar neutral hydrogen and dust attenuation towards SN 
1987A. Scaling the Zhekov et ah (2006) model X-ray spectrum to the observed day 9885 
0.3 - 8.0 keV flux, the total 0.01 - 8.0 keV luminosity is Lx (0.01 - 8.0 keV) = 5.6 x 10^® erg 
s“^. Assuming that the only source of EUV/soft X-rays is that observed by Chandra and 
that the neutral hydrogen in the outer ejecta intercepts half of the soft X-rays emitted by the 
ring, the Lyman-a heating efficiency must be ~ 60% to reproduce the F122M image of SN 
1987A. However, slower shocks that contribute more to the optical/UV emission of the ER 
than to the X-rays will also be a source of EUV irradiance of the outer debris (Fransson et ah 
2013), meaning that the actual soft X-ray heating efficiency is likely substantially lower. 

The HST observations presented here were acquired as part of the Cycle SAINTS 
program (13401 and 13405). 
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Fig. 1. — 20 August 2014 STIS G750L spectra with WFC3 filter curves: dashed cyan lines 
represent F645N and F665N, while the solid red line represents the F656N hlter. The F658N 
hlter is shifted about 2 nm to the red of F656N and is not shown here. Narrow emission 
lines from the ER are labeled in yellow and the broad RS arcs are labeled in orange. The 
dotted purple line shows the position of the equatorial ring. These spectra have an effective 
slit width of 0.5" in the East- West direction of the images shown in Figure 2-4. 
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Fig. 2. — The left column shows the original F645N (“H-a Blue”) and F665N (“H-a Red”) 
images (21 June 2014), with the blow up region outlined in gray. The right column shows an 
enlarged version after subtracting emission from the equatorial ring. The [Nil] contours are 
shown in magenta to identify the location of the ER. The regions of interest are approximately 
outlined in dashed white lines. 
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Fig. 3. — (top left) 2011 COS spectrum, with F122M filter curve overplotted (dashed red 
line). At upper right, R{La/Ha) (photon flux ratio) in a 0.225" E-W spatial extraction 
region (orientation indicated on image below) with normalized X-ray flux in blue circles. 
Emission from the ER prevents a measurement of R{La/Ha) at angles > ± 0.6". R{La/Ha) 
is elevated towards the regions of brightest X-ray emission. The ACS/SBC F122M Lyman-a 
image is shown at bottom (June 21 2014). 
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Fig. 4. — ACS/SBC F122M Lyman-« image with 19 March 2014 Chandra (0.3 - 8.0 keV) 
contours overlaid in blue (top). Below, spatial prohles of the western edge of the ER region 
show that the Lyman-a emission arises interior to the equatorial ring as traced by X-ray 
(blue), [O III] (green), and low- velocity H-a emission (black). 
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diagonal lines connecting the filled circles represent the circnmstellar rings. The red (bine) 
vertical lines delineate the slabs in which H-a emission from the reverse shock is visible 
throngh the “H-a Red” ( “H-a Bine”) filter. The A, B, C, D labels refer to the image 
locations in Fignre 2; the “b” and “r” refer to the “H-a Bine” and “H-a Red” images, 
respectively. Regions colored bine and red appear in Fignre 2, while regions colored gray are 
invisible to the H-a imaging. 


